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We determine the interface tension of the finite-temperature electroweak phase tran- 
sition in a numerical investigation of the SU(2)-Higgs model on a four-dimensional 
lattice with temporal extension Lt = 3. In this simulation the chosen parameters cor- 
respond to a Higgs boson mass of about 16 GeV. As a result the interface tension shows 
only small scaling violations in comparison with previous studies for Lf = 2 lattices. 



We also report on some experiences with autocorrelations in the applied Monte Carlo 
simulations of two-phase systems. 

1 Introduction 

According to the ideas of cosmology and elementary particle physics the electroweak phase 
transition takes place with decreasing temperature between a symmetric phase and a phase 
with broken symmetry, the Higgs phase. An order parameter of this transition is the vacuum 
expectation value of the scalar field, which has a non-zero value in the broken phase and 
vanishes in the high-temperature, symmetry- restored phase 0. 

The main interest in the electroweak phase transition in the early universe emerges from 
the question, whether it alone can provide a mechanism for the observed baryon asymmetry 
within the minimal Standard Model |§. A necessary condition for this scenario, essentially 
depending on the Higgs boson mass value tuh, is a transition of strong enough first order 
type, whereas for a weak transition every asymmetry generated at the phase transition would 
be washed out in the Higgs phase, because at temperatures T larger than the vector boson 
mass my/ the baryon and lepton number violation in the Standard Model is enhanced. 
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After the neglection of the U(l)-gauge field and the fermionic sector one is left with 
the SU(2)-Higgs model containing all important ingredients for these phenomena. The 
common theoretical tool for its examination is resummed perturbation theory [|], |5], |6j , which 
works well for low and intermediate Higgs masses (m# < 50 GeV). Besides the known 
problems caused by the infrared divergencies in the symmetric phase, its validity is not 
ensured any longer for higher values of where the strength of the electroweak phase 
transition is expected to decrease rapidly. This motivates a nonperturbative treatment on 
space-time lattices in four dimensions as well as in reduced three-dimensional models, see 
e.g. refs. |7], 14, 16, 17] and ref. for a more general review of the whole subject. 



The present work is a completion of numerical simulations of the four- dimensional SU(2)- 
Higgs model on the lattice in ||, [HJ, which focussed on Higgs boson masses below 50 GeV. 
Since this mass range is ruled out by the actual experimental bound of m H > 65 GeV, 
their main goal is to establish suitable methods for extracting physical quantities from the 
lattice investigations, and to hint at systematic errors in resummed perturbation theory or 
dimensional reduction approaches. 

In this publication we deal with the interface tension a. This quantity plays a prominent 
role in the course of the electroweak phase transition, because its magnitude is a measure 
for the strength of this transition. The nucleation rate of the Higgs phase in the symmetric 
phase is given in the thin-wall approximation by |2| 

Here T c denotes the critical temperature, Ae the latent heat, and r\ = is the so-called 
supercooling parameter. The prefactor r can be approximated by T 4 ||. If cx 3 /(Ae) 2 T c is 
large, signalling a strong first order phase transition, a substantial supercooling has to be 
expected, which would lead to an additional suppression of the sphaleron rate. 

In the following we use a low value of the Higgs boson mass, where the phase transition 
is quite strong. For the determination of a the two-coupling method in the scalar hopping 
parameter is employed. Compared to transfer-matrix techniques and histogram methods, it 
gives an optimal ratio between desired accuracy and required CPU-time for the SU(2)-Higgs 
model 00,0. 

2 Lattice simulation 

The lattice action of the SU(2)-Higgs model is conventionally [Tj| parametrized as 

S[UM = /^(l-lTWp) 
p 

+ E { l Tr feW + A [|Tr (<p+<p u ) - l] 2 - K E Tr {<pUfP**<P*) } ( 2 ) 
xen ^ ' v=i ' 
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in terms of the SU(2)-link variables U Xjfi and the site variables (p x , which are com- 
plex 2 (g) 2 matrices, representing the gauge and scalar degrees of freedom, respectively. 
Up = U x ^U x+ ^ u lI^ +0 JJ XjV is an elementary plaquette, and we often decompose the Higgs 
field as ip x = p x a x with p x G R >0 and a x G SU(2). If not stated otherwise, the lattice con- 
stant a is assumed to be a = 1, and the lattice volume is denoted by Q. The identifications 
g 2 = 4//3, wig = (1 — 2X)/k — 8 and Ao = A/4k 2 relate the lattice parameters (3, k and A to 
the bare gauge coupling, scalar mass and quartic coupling of the corresponding continuum 
theory. 

A suitable observable for the interface tension is the density of the yj-link operator L^^, 
which is conjugate — in the thermodynamic sense — to the hopping parameter k and is 
itself an order parameter of the phase transition: 

! 4 

Lcp ^ 4Q ^ ^ Lvvxtl ' Lvvxtl ^ i^+f^ Ux ^ x ) ■ ( 3 ) 

When simulating finite-temperature field theory, one utilizes lattices with spacelike ex- 
tensions much larger than the temporal extension: L s ^> L t , s G {x, y, z}. The physical 
temperature is given by the timelike lattice extension via T = l/aL t , and the approach to 
the scaling region of the model in the continuum limit is realized as a —>■ with T fixed, 
i.e. as L t — > oo. This limit goes along the lines of constant physics, on which renormalized 
couplings and masses are held fixed and only the lattice constant a is varying. 

In this spirit we extend the measurement of a at a -1 = 2T C in ref. |10[ to a -1 = 3T C with 
the purpose of gaining control over possible lattice artifacts. The parameters in the former 
analysis were (3 = 8.0 and A = 0.0001, leading to a renormalized gauge coupling of g 2 R ~ 0.56 
and a Higgs mass of m# ~ 16 GeV. The physical mass scale is set by the vector boson mass 
value raw = 80 GeV at T = in lattice units. A two-coupling simulation on a lattice of size 
2 x 16 2 x 128 resulted in a phase transition point at k c = 0.12830(5) and a finite-volume 
estimator for the interface tension of <r/T c 3 = 0.84(16). 

When passing over to smaller lattice spacings, we have to scale all lattice extensions 
accordingly in order to keep the physical volume constant. In the case of low Higgs mass 
this is possible with an acceptable demand of computer resources. With increasing m# the 
situation becomes worse, because the phase transition weakens and thus one needs larger 
physical volumes to obtain a stable two-phase situation. Therefore we choose a lattice of 
size 3 x 24 2 x 192 with the changed parameters (3 = 8.15 and A = 0.00011. They have been 
obtained by an integration of the one-loop perturbative renormalization group equations with 



the transition point of the L t — 2 lattice as a starting value fpLOf] . As shown in ref. |fL2|| , a 
sufficiently precise estimate for the critical point in k c is only available by numerical methods. 

In our Monte Carlo simulations we use an optimized combination of heatbath and over- 
relaxation algorithms. In particular, the introduction of the simultaneous overrelaxation 
of all four cartesian components of the Higgs field [0, instead of the Higgs field length- 
overrelaxation as proposed in M, has substantially reduced the integrated autocorrelation 
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time Tint- In this way we gained factors of 3 — 10 compared to the values previously found 
IT] . This r int -behaviour is reflected in table [1| and figure [1| for the case of a two-phase 



in 



simulation, where as typical examples the autocorrelation functions T(t) of the <^-link oper- 
ators (§) in both phases and their difference AL V are considered. The observation that this 
difference has a significantly lower r int will become relevant in the next section. 



heatbath 


overrelaxation 


Tint in sweeps 




<P>x 






Px 


'■Px 


r(l) 
J-/tp 


r(2) 
Liip 




1 


4 


3 


3 


1 




20(6) 


9.4(3.3) 


9.2(3.1) 


1 


1 


1 






3 


6.3(1.1) 


2.2(4) 


2.0(4) 



Table 1: Autocorrelation times for a 2—k simulation. Each 
updating sweep consists of a sequence of different algorithms 
as given by the numbers in the left part of the table. 




t in sweeps t in sweeps 



Figure 1: Normalized autocorrelation functions for a 2—K simulation with p- 
overrelaxation (left) and (p-overrelaxation (right) in the updating sequence. 



3 Two-coupling method 

The calculation of the phase transition point and the interface tension is done by a modified 



version of the two-coupling method [19]. Since this method has been used in the SU(2)-Higgs 

ylj , we only sketch its main idea here. One takes a periodic 



model before, see refs. 
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lattice with one elongated spacelike direction, i.e. L z ^> L x = L y ^> L t , and divides the 
corresponding lattice volume into two halves with different scalar hopping parameters 

k = (ki, k 2 ) = < k c for z < L z /2 , k 2 > k c for z > L z /2) . (4) 

Hence the lower half is forced in the symmetric phase and the upper one in the Higgs 
phase. If in the foregoing inequality the transverse directions L x and L y are sufficiently 
large, the system resides in a mixed-phases state and gives rise to an interface pair at the 
phase boundary perpendicular to the z-direction. 

In a first step one initializes a two-phase situation by a simulation with ^-values far 
away from the transition point, whose location can roughly be determined by hysteresis 
runs. Subsequently, the distance in k is more and more diminished, and the smallest k- 
interval, for which the system still resists to turn over into one single phase, gives lower and 
upper bounds for the critical k. On the 3 x 24 2 x 192 lattice we obtain the estimate 

K e = 0.128110(3) . (5) 

The interface tension is defined as the free energy F per unit area of the walls separating 
the two phases. As derived in ref. TO], its lattice version 



a 3 cr 



^— {f( Ki , k 2 ) - \F{k u «i) - iF(« 2 , k 2 )\ (6) 



is related to the expectation values L^\kx, k 2 ) and L^\ki, k 2 ) of L v in each phase by 

a 3 a = - lim lim { («i - k 2 ) lim L z ■ AL v («i, k 2 ) 1 , (7) 

where AL v (ki, K 2 ) denotes their difference L^\ki, k 2 ) — L^\ki, k 2 ). With the (N + 2)- 
parametric Laurent ansatz 



N 

k 2 ) = 5— + V 7 ^(k, - + O (\k, - k c \ n+1 ) , i = 1, 2 . 



j=0 

this leads to the finite- volume estimator for the interface tension 

a 3 a = L z (c x + c 2 ) . (9) 

The inverse- linear term in eq. (H) is motivated by the fact that for Are = |/c$ — k c | 1 the 
free energy change between the two phases behaves as AF ~ O(Ak). The probability p for 
an interface being at zo > z is essentially given by exp{— const • Ak(z — zo)} and therefore 
J dzp ~ 1/Ak. 

Thus one proceeds in a similar way as for the /^-determination, but due to practical 
limitations on L z , one has to prevent the interfaces from touching each other by a choice 
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lattice: 3x24x24x192, fi = 8.15, A = 0.0001 1 



40 - 



20 



k,=0. 12776 



K,=0. 12846 



40 



20 



lattice: 3x24x24x' 92, (3 = 8. \5, A = 0. 00011 



/c, = 0. 12801 



Figure 2: Two-phase profiles of the z-slice expectation value L^z) of L v for our largest 
and smallest K-intervals. 



of large enough K-intervals. Figure |2| displays typical two-phase structures from the 2—k 
method. Note that the plateaus become narrower for smaller K-intervals, especially in the 
broken phase, but the phases are still clearly developed, and the interfaces continue to exist. 

The results of our simulations are presented in table |2|. Owing to correlations between 
L$ and \ the statistical errors on AL^ are usually smaller than those of . This was 
already suggested by the autocorrelations in figure [l] and is caused by shifts of the interfaces 
during the simulation, which cancel out in AL^ to some extent. 



Ki 


«2 


sweeps 


L (1) 
j-iip 


r(2) 


AL V 


0.12776 


0.12846 


5000 


1.5161(24) 


37.8216(70) 


36.3055(67) 


0.12781 


0.12841 


5000 


1.5149(26) 


34.7287(68) 


33.2138(60) 


0.12786 


0.12836 


5000 


1.5285(29) 


31.5085(77) 


29.9799(76) 


0.12791 


0.12831 


10000 


1.5398(19) 


28.1144(55) 


26.5746(57) 


0.12796 


0.12826 


10000 


1.5616(30) 


24.5264(76) 


22.9648(71) 


0.12801 


0.12821 


20000 


1.6019(32) 


20.5877(66) 


18.9858(59) 



Table 2: Results for L^', Ly and AL V on a 3 x 24 2 x 192 lattice. 
The errors are obtained by binning. 



In order to give a reliable estimate for the statistical error when fitting these 95-link aver- 
ages to the functions in eq. (§), we take their correlations into account by a bootstrap analysis 
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|20|j . The characteristic ingredient is to calculate secondary quantities from bootstrap sub- 
samples randomly taken with repetition from the original measurements, which itself can be 
interpreted as the empirical probability distribution of the observables under consideration. 
The errors of the fit parameters are extracted as half of the central 68.3%-interval of their 
distributions from fits performed on these sample averages. A more extensive description of 



this method in the same context is contained in ref. ITT 



We quote our final result for two four-parameter fits with x 2 -values of 2.58 and 0.82. 
Using (|D|), we find 

(10) 



a 



T 3 , 

c / L 



— = 0.764(52 + 47) 



whereby aT c — 1/L t = | in lattice units. The error consists of two parts, coming from a 
bootstrap analysis with 10000 iterations and the uncertainty of k c in eq. (HI). 



lattice: 3x24x24x192, /? = 3.15, A = 0.0001 1 




40.0 



lattice: 3x24x24x'92 ; fi = S.15, A = 0.0Q011 



0.12785 0.12795 
k<k =0.1281 1 



20.0 u 

0.1282 



0.1283 0.1284 
K>K„=0.1281 1 



Figure 3: Four-parameter least-squares fits of i = 1,2, separately in each phase 



At this place some explanations about the number of necessary fit parameters are in 



order. In the examinations of L t = 2 lattices ||T0| , |TTJ a three-parameter fit was sufficient 
to give a reliable value for the interface tension; an inclusion of higher order terms gave 
no improvement. The situation for L t = 3 seems to be different for the following reason: 
Basically, one has to keep in mind that the expansion in eq. (|^) is only a phenomenological 
ansatz. Furthermore the lattice volume Q has increased, and the phase transition, which 
in a strict sense is only present in infinite volumes, is more pronounced. Consequently, the 
contribution of the higher Laurent coefficients becomes more important for the slopes at 
larger K-intervals. 

In fact, the three-parameter fit has no satisfactory x 2 . A five-parameter fit gives 
ct/T 3 = 0.85(20 + 5) and x 2 -values equal to 1.42 and 0.82. This is compatible to ([10]) 



7 



within errors, but fairly sensitive to the number of fitted data points. A more careful inspec- 
tion of the bootstrap calculations reveals that the last fit parameter is not very significant in 
such cases, where its bootstrap error is roughly as large as its value; this holds true for the 
five-parameter fits and also the four-parameter fit in the symmetric phase. So for the sake 
of completeness we combined the three-parameter fit in the symmetric phase with the four- 
parameter fit in the broken phase to <5"/T c 3 = 0.793, although we are nevertheless convinced 
of the four-parameter fit in both phases to lead to the most reasonable result. 

The preceding remarks on the <£>-link correlations should have made clear that AL V is 
the most natural variable for estimating the interface tension via eq. (0). This requires that 
the chosen ^-intervals are symmetric with respect to k c . A four-parameter fit to an ansatz 
similar to fl8|) for AL V yields 



with x 2 — 0.39 and is illustrated in figure |4]. Here we only quote the statistical error, which 
now comes from 1000 normally distributed random data, since the different AL^-averages 
are uncorrelated. Note the perfect agreement of this result, and especially of its error, with 
the numbers in eq. ( |T0"D above. 




(11) 



attice: 3x24x24x192, (S = 8A5, A = 0.0001 1 



35.0 



25.0 



20.0 



4 



0.0001 5 



0.00035 



0.00055 



0.00075 



Figure 4: Four-parameter least-squares fit of AL^ as a 
function of — K\ . 
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4 Discussion 



We have determined the interface tension of the four- dimensional SU(2)-Higgs model with 
the two-coupling method. All estimates from acceptable fits with a reasonable number of 
parameters show a very good consistency. 

To our knowledge this is the first attempt to supply information about the scaling of o 
when going over to a finer lattice. Compared to the L t = 2 results 

4) =0.84(16), (^) =0.83(4) (12) 

1 c J 2-K V c / hist 

from refs. || |10[ — the second number referring to the histogram method — the accuracy of 
the 2-k estimate has been improved to the value cx/T c 3 = 0.76(10) in eq. (0). The observed 
small deviation between L t = 2 and L t = 3 confirms the smallness of scaling violations, as 



e.g. also found recently for the critical temperature in |12 



Finally, we confront our result with the perturbative estimate || up to order # 4 ,Aq 

£s) =0.78(1), (13) 

c / pert 

with an error coming from the uncertainties in the renormalized parameters on the lattice. 
As expected, the agreement with perturbation theory on a quantitative level in this range 
of the Higgs boson mass is excellent. 
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